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In this paper, current-induced spin polarization for two- 
dimensional electron gas with a general spin-orbit in- 
teraction is investigated. For isotropic energy spectrum, 
the in-plane current-induced spin polarization is found 
to be dependent on the electron density for non-linear 
spin-orbit interaction and increases with the increment 
of sheet density, in contrast to the case for fc-linear 
spin-orbit coupling model. The numerical evaluation is 
performed for InAs/InSb heterojunction with spin-orbit 
coupling of both linear and cubic spin-orbit coupling 
types. For J-type short-range electron-impurity scatteri- 



1 introduction Spintronics, where the spin degree 
of freedom is manipulated to control the electronic de- 
vices, has been becoming a rapid field of condensed matter 
physics 0]. Current-induced spin polarization (CISP) dis- 
closes the possibility of the spin polarization generated 
in semiconductor directly by the electric field. It refers to 
a spatially homogeneous spin polarization in two dimen- 
sional systems due to an in-plane charge current. CISP 
due to spin-dependent scattering was first reported by 
D'yakonov and Perel' in 1971 J2]. Later it was realized 
that such transport phenomenon exists in semiconductor- 
based two-dimensional electron gas (2DEG) without struc- 
ture or bulk inversion symmetry 13)4.(5). Experimentally, 
CISP was first measured by Silov et al. in two dimensional 
hole system with the help of polarized photoluminescence 
technology [6j. Observations of CISP in strained semi- 
conductors have been reported by Kato EM). And later 
Sih et al. demonstrated the existence of CISP in AlGaAs 
quantum well J9) . 

So far, the theoretical investigations about the CISP are 
mainly focused on the 2DEG with fc-linear SOC, such as 



ng, it is found that the current-induced spin polariza- 
tion increases with increasing the density when cubic 
spin-orbit couplings are considered. However, for remote 
disorders, a rapid enhancement of current-induced spin 
polarization is always observed at high electron den- 
sity, even in the case without cubic spin-orbit coupling. 
This result demonstrates the collision-related feature of 
current-induced spin polarization. The effects of differ- 
ent high order spin-orbit couplings on spin polarization 
can be comparable. 
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Rashba SOC due to structure inversion asymmetry Hr = 
a(k y <j x — k x (jy ) ll3ll5l[T0l, linear Dresselhaus SOC due 

to bulk inversion asymmetry Hjy = (3{k x a x — k y a y ) 
lUD . and the combination of linear Rashba and Dressel- 
haus coupling types |12lll3lfT4l . Here er = (<r x , a y , a z ) 
represents the set of Pauli matrices, and fc = (k x , k y ) = 
fc(cos#, sin(9) is the two-dimensional momentum, a and 
(3 are linear Rashba and Dresselhaus SOC factors, respec- 
tively. CISP is found to be proportional to the SOC con- 
stant and independent of the electron density for 2DEG 
with fc-linear Rashba or Dresselhaus SOC when a dc elec- 
tric field is applied l3ll5lfTTl. 



However situations may become different for 2DEG 
with non-linear SOC, just like the spin Hall effect lfl5l 
fTMT7lfTMT9ll20ll2T1l2a. the investigation of CISP on this 
system is desirable. Liu et al. studied the CISP in hole- 
doped two dimensional system lacking structure inversion 
symmetry, and they found that CISP is dependent on the 
Fermi energy, in vivid contrast against the case in 2DEG 
with fc-linear SOC [23). While it was pointed by Liu et 
al. that the "spin" for hole system is actually the total an- 
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gular momentum, where the spin of hole system is not a 
conserved physical quantity fl23l . Hence, we expect the 
behavior of the spin polarization of conduction electron 
in the presence of non-linear SOC. We know that at high 
electron sheet density, the cubic term of Dresselhaus SOC, 

Hp = r}(k m ky<J x — k x k y a y ) 11241 . (77 is the cubic Dres- 
selhaus SOC constant), has to be taken into account l25l . 
Recently, by using the double-group representations, Car- 

toixa et al. Il26l found that there is another cubic term for 

(3) 

heterojunction due to bulk inversion asymmetry, H y B f A = 
C{kyPy — k x a x ), (i. e. the last term of Eq. (3) in Ref. 
1261 ). Here C = a 2 (3/6 is the cubic SOC constant with 
a as the well width. Apart from the type of SOC, most 
of the pioneering researches on CISP treat the electron- 
impurity collision with a simple momentum-independent 
form, described by a relaxation time r. However in realis- 
tic 2DEG, the electron density is not large enough to screen 
the charged impurities. The interaction between electron 
and disorder is long ranged. 

In this paper we consider the CISP in 2DEG with a 
general SOC, which can be applied to describe Rashba, 
linear and cubic Dresselhaus SOCs, and many other SOCs. 
For the simple isotropic energy band form, the analytical 
result of CISP is obtained. The numerical evaluation for 
the electron system with both linear and cubic SOCs due 
to bulk inversion asymmetry is also performed, considering 
both short- and long-range disorders. 

2 formalism We consider a two-dimensional non- 
interacting electron system with a general SOC, described 
by the following one-particle Hamiltonian: 



Dresselhaus SOC, a = 0, A = 0, Mi = 1, M 3 = -1; and 
the case & = -\ir), A = 1, M x = 3, M 2 = 2, M 3 = 1 
corresponds the cubic Dresselhaus term. 
With the help of the local unitary matrix 



Ut 



V2 (1 



1 



where \k satisfies 



tanxfe = 



a, sin M 3 9 — a r cos M 3 9 
014 cos M36 + a r sin M36 ' 



(3) 



(4) 



the Hamiltonian (fTJ can be diagonalized into H = 
diag[ei(fc), s 2 (k)] in the helicity basis. Here 



e„(fc) = e (fc) + (fc), 



(5) 



with e M {k) = \a\k Ml (sin M 2 9) x and p, = 1, 2 as the he- 
lix band index. We note that when A = 0, the energy dis- 
persion e^(fe) becomes isotropic function of wave vector 
k, whereas for the case A = 1, the energy spectrum relies 
on the angle of momentum. 

When the system is driven by a weak dc electric field 
applied along the x direction, E = Eq£. Following the 
procedure of Ref. l27l . the kinetic equation of the distri- 
bution function p(k) can be derived, with the equilibrium 
distribution function 



» P [ei(fc)] 
n P [e a (fc)] 



(6) 



H = e Q (k) + b(k) ■ er. 



(1) 



For simplicity we have assumed that £o(k), the energy dis- 
persion in the absence of SOC, is isotropic function of mo- 
mentum k. We set h = 1 throughout this paper. As a fur- 
ther simplification of SOC, we shall later specialize to the 
model, where the spin orbit field b(k) has the form 



b x {k) + ib y (k) = ak Ml (sin M 2 9) 



\\ e iM 3 6 



(2) 



Here the complex number a = a r + ion is the general cou- 
pling constant, irrespective of the momentum k, with a r 
and <Xi as the real and imaginary part, respectively. A = 0, 1 
and Mi, M 2 , M3 are integer numbers. It will be noted later 
that the index A determines whether the energy spectrum 
is isotropic. The number Mi usually is positive. It is found 
that when A = and Mi = M3, our model becomes the 
one in Ref. ||2T1| . where the model is used to discuss the spin 
Hall effect. Due to time reversal symmetry, the spin orbit 
field b(k) satisfies b(k) = —b(—k). Therefore for the case 
A = 0, ;. e. isotropic energy spectrum, M3 must be an odd 
integer number; while for the case A = 1, ;. e. anisotropic 
energy spectrum, M 2 + Ms must be an odd integer num- 
ber. Now we list some special SOC forms: for pure Rashba 
SOC, a = -ia, A = 0, Mi = M 3 = 1; while for fe-linear 



Here n-p (x) is the Fermi distribution. The distribution func- 
tion p(k) to first order of electric field comprises two 
terms. The first term is written as 



o (1) - o {1) - 



eE dxk 
'4e M (fe) dk x 



|n F [ei(fc)] - np[e 2 (fe)]|. 



(7) 

And the second term (k) is determined by the set of 
equations with the form 



eE 



dn F [£^(k)] 



dk x 



•E 



\u(k 



<7)| 2 ^' 



x [p%(k)~P%(q)] S(e,(k) 



4e M (fc)Rcpf 2 ) (fc) = 7T W fc - Q)\ 2 ^ 

qfifi' 



P { Zl{k)-p%,{q) S(e,(k) 



V (<?))> 
(8) 



(9) 



Here = 1 + (-1)^ cos(xfc - X q ) and = 



(— 1) M sin(xfc — x q ). Kepi 2 '(k) represents the real part 
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of the off-diagonal distribution function p^{k). u(k) is 
the electron-impurity scattering matrix. Note that the above 
equations of distribution function, Eqs. (HJ and (O, have 
been derived in Refs. ll3l4l[T3l. 

Finally, in helix spin basis, the single-particle opera- 
tors of spin polarization are given by Si = U^iaiUk with 
i = x,y, z. The corresponding macroscopical quantities 
are obtained by taking the statistical average over them, 
Si = J2 k Tr[p(k)Si], and are expressed as 



S x = ^^sinxfe[p 2 2(fe) - Pn(fe)], 



kft 



Sy = \ ^COS X k[p22{k) -p U (k)], 



ff, = 5^Repia(fc). 



(10) 

(11) 
(12) 



3 spin polarization 

3.1 analytical result For this angle-dependent SOC, 
the formulas to be derived will become much less trans- 
parent, and the integrals are more difficult to solve ana- 
lytically. Therefore, we first limit ourselves to isotropic 
energy spectrum, i. e. A = and the parabolic case 
£o(fc) = tt— ■ The spin polarization is examined in the 
presence of electron-impurity scattering with 5-potential, 
\u(k — q)\ 2 = riiUQ. Here m is the effective mass of 
two-dimensional electrons and n, is the impurity density. 

Keeping only the lowest-order of spin-orbit interaction, 
the diagonal elements of distribution function can be ob- 
tained analytically from Eq. ([S]). For M3 = ±1, the diago- 
nal elements of p^ (k) take the form 



eE Q T 



k + m\a\(l - M 1 )(2ttN) — 



x cosQ5(e k i - e F ), (13) 
eE r 



k-m\a\(l - M 1 )(2ttN)- 



cos68(sk2 - ef)- 



(14) 



And for the case M3 > 1 or M3 < — 1, they are given by 



pW(k) = -eE T 



dEi(k) 



-eE r 



dk 
de 2 {k) 



Ok 



cos95(s k i - E F ), (15) 
cos9S(e k2 - e f ), (16) 



with r = 1 /mriiUQ as the relaxation time, Ep as the Fermi 
energy. With the help of Eq. (O, the off-diagonal distribu- 
tion function can be obtained directly. Finally, we find the 



spin polarization takes the form 



Si, — 




(27riV) J 



|M 3 | = 1 
|M 3 | > 1 
\M 3 \ = 1 
|M 3 | > 1 



(17) 

(18) 
(19) 



We come to the conclusion that the in-plane CISPs exist 
only when the winding number |A/ 3 | = 1, nevertheless 
the out-of-plane component of CISP is always zero for 
this 2DEG with isotropic general SOC. As expected, for 
M 3 = ±1 and Mi = 1, we obtain the CISP for 2DEG 
with Rashba or fc-linear Dresselhaus SOC, in agreement 
with previous theoretical studies [3 5l lTT1fT2lfT3l . Our re- 
sults imply that, in contrast to the 2DEG with Rashba or k- 
linear Dresselhaus SOC, the in-plane CISPs for non-linear 
SOC system, Mi > 1, become dependent on the electron 
density and enhance for high density semiconductor. The 
above analytical calculation is valid in the weak SOC case, 
EM(kp) <C £o(A;f), but the relationship between £m(&f) 
and t _1 is arbitrary. Here kp is the Fermi wave vector. For 
strong SOC case, the dependence of CISP on the electron 
density can be evaluated numerically. 

It should be noted that the in-plane CISP comes from 
the interband processes, arising from the SOC, which can 
be seen from Eqs. ( [Tol l and (fTTl i. In the absence of spin- 
orbit interaction, pn = P22 leads to the vanishing CISP. 
Hence, although the result, we obtained here, is for the 
parabolic case, the nonparabolic contribution of Eo(k) to 
CISP may only change its value slightly through the Fermi 
energy. Further, it can be confirmed below by the numeri- 
cal calculation. In general, SOC field will be the combina- 
tion of Eq. <|2J with both linear and high order terms. From 
our analytical result, one can deduce that with increasing 
the high order SOC constant, the density-related feature of 
CISP will become more and more evident. 

3.2 numerical result Now we perform the numerical 
evaluation for CISP in InAs/InSb heteroj unctions without 
the additional large bias voltage, where the main SOC con- 
tribution terms arise owing to the absence of the center of 
inversion in the bulk material. 

Further, to take account of the nonparabolicity of the 
energy band of InAs, we use the isotropic Kane band 
model: 



e (k) = 




(20) 



where 7 w l/e s is the nonparabolic parameter, with e g as 
the energy gap between the conduction and valence bands. 
Note that Kane energy band becomes the parabolic case 
for vanishing 7. The nonparabolic factor in the numerical 
calculation is set to be 7 = 2.73 eV _1 for InAs |[28l . The 
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Figure 1 rnS x /Eo as functions of electron density for 
(a) short- or (b) long-range impurity scattering for 7 = 
2.73 eV" 1 . Here r) = 1.0 x lCT 28 eVm 3 , g is the effec- 
tive g-factor, and /is is the Bohr magneton. The thin solid 
line in (a) is obtained for 7 = and 77 = 770- The thin and 
thick solid lines in inset of (b) is calculated when rj = 770 
for 7 = and 7 = 2.73 cV -1 , respectively. The unit of 
electron density TV is 10 11 cm" 2 , and the unit of riiS x /Eo 
is 10 26 g^Es/Vm 3 . 



(5-form short-range or the remote charged impurity scat- 
tering is considered in the calculation. The scattering ma- 
trix of remote electron-impurity scattering takes the form: 
\u(q)\ 2 ~ n. l e~ 2sq I(q) 2 J29), with I(q) as the form fac- 
tor. We set the electron effective mass at the band bottom 
m = 0.04m e (m e is the free electron mass), remote impu- 
rities in InSb barrier are located at a distance of s = 10 nm 
from the interface of the heterojuction II22II28I . 

3.2.1 Hp^+Hp^ First we consider the system with 
linear and cubic Dresselhaus SOC. In this case, the spin 
orbit field b(k) = (/3k x + r)k x k 2 ,—fik y — rjk y k 2 ). At 
the same time, the equations about distribution function, 
from Eq. (O to Eq. (O, can be obtained, by substituting 
the new form of energy e&f(fe) and \k with £m(&) = 
y/b x (k) 2 + b y {k) 2 , Xk = ~ tan" 1 It is noted that 

the energy spectrum becomes anisotropic completely. 

From Eqs. ([8]) and ( [Tol l, we find that the x compo- 
nent of CISP is inverse proportional to the impurity den- 
sity m. Therefore, riiS x /Eo is plotted in these figures. In 
the calculation, we set linear Dresselhaus SOC coefficient 
(3 = 1.0 x 10~ n eVm. The short-range impurity scattering 
is considered with relaxation time r = 10 ps. When the 
cubic Dresselhaus SOC is considered, it is found that CISP 
is still along the x direction. 



Figure 2 Dependencies of riiS x / E$ on cubic Dresselhaus 
SOC parameter for (a) short- or (b) long-range impurity 
scattering. 



In Fig.Q] the x component of CISP is plotted as func- 
tion of electron density. For short-range disorder, the den- 
sity dependence of CISP can be observed when r\ 7^ 0. 
Note that the resultant S x is proportional to relaxation time 
r. With the increment of the sheet density, CISP increases 
monotonously, and almost saturates at high density for 
large rj. From Fig. [TJb), for long-range electron-impurity 
scattering it is evident that, unlike the case for short-range 
disorder, here CISP always increases with ascending the 
density even for the system without cubic SOC. In the 
parameter regime, N < 10 12 cm~ 2 , long-range disorders 
have strong effect on CISP, where CISP increases rapidly 
with the rise of density. It can be seen that the role of cu- 
bic term of Dresselhaus SOC on CISP becomes important 
at high sheet density for both short- and long-range colli- 
sion. CISP for parabolic energy band is also plotted in this 
figure with a thin line. We find that the weak effect of non- 
parabolicity on CISP appears at high density. 

CISP is shown as a function of the cubic Dresselhaus 
SOC parameter rj in Fig. |2] For momentum independent 
potential, CISP begins with the value em(3TEo/2ir, inde- 
pendent of the density, and increases with ascending r\. In 
Fig. EJb), the calculated CISP for long-range collision is 
almost linear proportional to cubic Dresselhaus constant rj. 

3.2.2 + + H^] A In this subsection, the 

additional high-order contribution H B j A due to bulk in- 
version asymmetry in Ref. l26l is also considered. Now 
the spin orbit field becomes b(k) = ((3k x + rjk x k 2 — 
(k^,, — [iky — ijkyk 2 + (ky), and the corresponding £jv/(fe), 
Xk can be obtained analogously. We take the well width 
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N ( 10" cm 2 ) 

Figure 3 riiS x /Eo is shown as functions of electron den- 
sity for (a) short- or (b) long-range impurity scattering 

(3) 

when the additional cubic SOC term H B } A is considered. 
The thin solid lines in (a) and (b) are obtained when £ = 
and 77 = jyo for short-range and long-range collisions, 
respectively. The inset in (b) shows the dependencies of 
riiSx/ 'Eq on N at high density regime. The other parame- 
ters are the same as those in Figure [TJ 



a = 5nm. The calculated riiS x /Eo as functions of elec- 
tron density N is shown in Fig. [3] 

(3) 

When the additional high-order term H B } A is in- 
cluded, the magnitude of CISP rises for both short- and 
long-range disorders. For short-range scattering, riiS x /Eo 
always increases with ascending the density. However, at 
high density, the magnitude of CISP for large 77 may be 
less than the one for small 77. This is due to the interplay 
between two cubic terms. It has been seen that CISP satu- 
rates at high density when only the Dresselhaus cubic term 

(3) 

H D is included. However, one can see, from the thick 
solid line in Fig. Oa), this behavior will not occur when 

(3) 

we only include the term H£} A . For large 77, the effect on 

(3) 

CISP of cubic Dresselhaus term H D exceeds the one of 

(3) 

this additional cubic term H B j A . CISP saturates again at 
high density, hence its magnitude becomes less than the 
one with small rj. However, such phenomenon can not be 
observed for the case of long-range collision. 



4 conclusion In summary, the CISP for 2DEG with 
a general SOC is investigated. For isotropic energy band, 
we find that the in-plane CISP becomes density-dependent 
for non-linear SOC, and increases with enhancing the sheet 
density. We have numerically studied the linear and cubic 
SOC contributions to CISP, considering both the short- and 
long-range disorders. For short-range collision, we have 
demonstrated the dependencies of CISP on density when 
high-order SOCs are included. When impurity scattering 
becomes long ranged, however, CISP increases rapidly 
with raising the density even for the system without cubic 
SOC. Our investigation indicates that the remote disorder 
has a strong influence on spin polarization, and the magni- 
tude of CISP strongly relies on the scattering matrix. The 
contributions of different cubic SOCs to CISP can be com- 
parable. 
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